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Abstract: In the present work, we study and analyze an efficient iterative coupling method for 
a dimensionally heterogeneous problem . We consider the case of 2-D Laplace equation with non 
symmetric boundary conditions with a corresponding 1-D Laplace equation. We will first show 
how to obtain the 1-D model from the 2-D one by integration along one direction, by analogy with 
the link between shallow water equations and the Navier-Stokes system. Then, we will focus on 
the design of an Schwarz-like iterative coupling method. We will discuss the choice of boundary 
conditions at coupling interfaces. We will prove the convergence of such algorithms and give some 
theoretical results related to the choice of the location of the coupling interface, and the control of 
the difference between a global 2-D reference solution and the 2-D coupled one. These theoretical 
results will be illustrated numerically. 

Key-words: dimensionally heterogeneous coupling, domain decomposition methods, multiscale 
analysis 
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Une methode de couplage de type Schwarz dans le cadre 
d'un probleme multi-dimensionnel 

Resume : Dans ce document nous etudions et analysons et une methode de couplage multi- 
dimensionnel iterative. Nous considerons le cas de l'equation de Laplace 2-D avec des conditions 
aux bords non symetriques, couplee avec une equation de Laplace 1-D correspondante. dans 
un premier temps nous montrons comment obtenir le modele 1-D a partir du modele 2-D par 
integration verticale et par analogie avec la derivation des equations de Saint- Venant a partir 
des equations de Navier-Stokes. Ensuite nous presentons un algorithme de couplage de type 
Schwarz. Nous discutons le choix des conditions aux interfaces de couplage. Nous demontrons 
la convergence de tels algorithmes et donnons quelques resultats theoriques sur le choix de la 
position des interfaces de couplage. Un resultat theorique sur le controle de l'erreur entre la 
solution globale 2-D de reference et la solution 2-D couplee sera aussi donne. Enfin nous illustrons 
ces resultats numeriquement. 

Mots-cles : couplage multi-dimensionnel, decomposition de domaine, analyse muli-echelles 
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1 Introduction 

Hydrodynamical phenomena can be described by a wide variety of mathematical and numerical 
models, spanning a large range of possible levels of complexity and realism. When dealing with 
the representation of a complex fluid system, such as an ensemble of rivers and channels or a 
human blood system, the dynamical behavior of the flow is often spatially heterogeneous. This 
means that it is generally not necessary to use the most complex model everywhere, but that 
one can adapt the choice of the model to the local dynamics. One has then to couple several 
different models, corresponding to different areas. Such an approach is generally efficient from 
a computational point of view, since it avoids heavy computations with a full complex model in 
areas where a simpler model is able to represent the dynamics quite accurately. Thus this makes 
it possible to build a hybrid numerical representation of an entire complex system, while its 
simulation with a unique model would be either non relevant with a simple one or too expensive 
with a complex one. 

In such a hierarchy of models, the simplest ones are often simplifications of the more complex 
ones. Let mention for instance the so called "primitive equations", which are widely used to 
represent the large scale ocean circulation, and are obtained by making some assumptions in 
the Navier-Stokes equations. It is important to note that such simplifications may involve a 
change in the geometry and in the dimension of the physical domain, thus leading to simplified 
models which are m-D while the original one was n-D, with n > m. An obvious example is 
given by the shallow water equations, which are derived from the Navier-Stokes equations by 
integration along the vertical axis, see for instance [T] for a rigorous mathematical derivation using 
asymptotic analysis techniques in the 2-D to 1-D case . Such a coupling between dimensionally 
heterogeneous models has been applied for several applications. Formaggia, Gerbeau, Nobile 
and Quarteroni [5] have coupled 1-D and 3-D Navier-Stokes equations for studying blood flows 
in compliant vessels. In the context of river dynamics, Miglio, Perotto and Saleri [3J, Marin and 
Monnier [3], Finaud-Guyot, Delenne, Guinot and Llovel [5], Malleron, Zaoui, Goutal and Morel 
[6] have coupled 1-D and 2-D shallow water models. Leiva, Blanco and Buscaglia [7] present also 
several such applications for Navier-Stokes equations. 

Several techniques can be used to couple different models, either based on variational, al- 
gebraic, or domain decomposition approaches. In the context of dimensionally heterogeneous 
models, in addition to the previously mentioned references, there exists also a number of papers 
on this subject for purely hyperbolic problems (e.g. [5J[H], [HI]), but we will not elaborate on 
these studies since our focus is much more on hyperbolic/parabolic problems. 

In our study, we will focus on the design of an efficient Schwarz-like iterative coupling method. 
The possibility of performing iterations between both models, i.e. of using a Schwarz method, is 
already considered in Miglio, Perotto and Saleri [3] and Malleron, Zaoui, Goutal and Morel [§]. 
This kind of algorithm has several practical advantages. In particular, it is simple to develop and 
operate, and it does not require heavy changes in the numerical codes to be coupled: each model 
can be run separately, the interaction between subdomains being ensured through boundary 
conditions only. These are important aspects in view of complex realistic applications. 

Our final objective is to design an efficient algorithm for the coupling of a 1-D/2-D shallow 
water model with a 2-D/3-D Navier-Stokes model. As a first step in this direction, the present 
study aims at identifying the main questions that we will have to face, as well as an adequate 
mathematical framework and possible ways to address these questions. We will perform this pre- 
liminary stage on a very simple testcase, coupling a 2-D Laplacian equation with a corresponding 
simplified 1-D equation. Seemingly similar testcases were addressed by Blanco, Discacciati and 
Quarteroni [TT] and Leiva, Blanco and Buscaglia Q~2] , but with different coupling methodologies 
(variational approach in |11| and Dirichlet-Neumann coupling in |12|). Moreover we have chosen 
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to use non symmetrical boundary conditions in our 2-D model, in order to develop a fully two 
dimensional solution, and our 1-D model is obtained by integration of the 2-D equation along 
one direction, by analogy with the link between the shallow water system and the Navier-Stokes 
system. The rigorous mathematical derivation of the 1-D model clearly highlights its validity 
conditions. 

Section 2 is devoted to the presentation of the 2-D Laplacian model, and to the derivation of 
the corresponding reduced 1-D model. Then a Schwarz iterative coupling algorithm is presented 
in Section 3, and its theoretical convergence properties are analyzed. In particular, the influence 
of the interface location is discussed. Finally numerical tests are presented in Section 4, which 
fully validate the previous analytical results. 



2 Derivation of the reduced model 

We are interested in the following boundary- valued problem in a domain f2 C M 2 : 



with a and k are nonnegative numbers that allows for Dirichlet, Neumann or Robin boundary 
conditions. 

In this section, we want to take advantage of the shallowness of the domain f2 (or some subdomain 
of fi) to derive a reduced model. As it is done for the derivation of the shallow water model (see 
e -9- OP)) we want to replace the complete 2-D model ([T]) by a (simpler) 1-D equation wherever it 
is possible (and keep the original 2-D model everywhere else) . We will finally obtain the coupled 
1-D / 2-D system ()15[) - (|16)) which we will analyze and simulate in the coming sections. 
In order to discriminate between 1-D and 2-D regions, we introduce the following definition: 

Definition 1. Let Clio be the subset of Q in which 2-D effects may be neglected, and f^D == 
f2\£7i£) the subset of Q in which 2-D effects cannot be neglected. 

Remark 1. Naturally the definition of flio depends on several features, such as the domain 
aspect ratio, the considered system of equations, forcing terms (including boudary conditions), 
etc. 

For the sake of clarity (see Figure [T]), we will assume that there exists H and L\ such that 



Let us now consider equation ([T]) with the following boundary conditions (see Figure [T] for 
the notations): 




(lb) 



(la) 



fl 1D = fin {x < Li} = (0,£i) x (0,H) and il 2D = On {a; > L x }. 



— Au(x, z) 



F(x,z), V(z,z)en, 



(2a) 




o, v(i,z)er T , 



(2b) 



< 



— (x, z) + ku(x, z) 

u(x, z) 
u(x, z) 



0, V(a;,z)er B , 

7i(.t,z), M{x,z) e r L , 

j 2 (x,z), V(x,z) e T R . 



(2c) 




(2e) 
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Figure 1: Typical computational domain £1, including a zone with a true 2-D behavior (x > L\) 
together with a shallow zone where we intend to use a 1-D model (x < L\). 



In order to derive the 1-D model in ilm = Q H {x < Li}, we introduce the following 
dimensionless variables and numbers: 



H 
I? 



(3) 



x z u(x, z) ~ L\ 

x=—, z=—, u[x,z) = — — — , F(x, z) = F(x, z) — and k = kLi, (4) 

L\ H U U 

where L\ (resp H) is the characteristic length (resp height) of Clio, £ is called the aspect ratio, 

and U is a characteristic value for u(x, z). 

The nondimensional form of equations ([2} in VL\d readtQ: 



d 2 u 1 d 2 u 

dx 2 e 2 dz 2 



F in VLid 



^ = o on ri = r T n dn 1D 

oz 

- -^ + ku = o on ro = r B n50i D 

£ OZ 



7i 

u 



on rl=T L . 



We assume (see [12] for the scaling of k) that 



(5a) 
(5b) 
(5c) 

(5d) 



F = 0(1), = Oil) and K = 0(e), 
ox 2 



(6) 



which is a sufficient condition to ensure that the 2-D effects are negligible in Hid- Indeed we 
deduce from equation (|5aj) that 



1 Since £1\d = (0, Li) X (0, H) we have n = ±e z in Il5bt and H5c> , 
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By vertical integration on (z, 1), and accounting for the boundary condition (|5b|) . we find: 

g=0<a (8) 

and finally: 

u(x,z) = u{x,0) + O(e 2 ). (9) 
Going back to original variables, we have: 

u(x, z) = u(x, 0) + O(e 2 ), V(x,z) e n 1D . (fO) 

We now introduce the averaging operator in the vertical direction. For any function / of z, 
we set: 

1 ' H 



f= Hj Q f[z)dZ - (U) 
We integrate equation §§§ for z G (0, iJ) and obtain: 

u(x) = w(x,0) + 0(e 2 ), Vx G [0,i]. (12) 

We now average equation ([5a} in the vertical direction, taking into account the Robin boundary 
condition (|5c[) on T l Bl and find: 

For every a; G [0,Li] we may use approximation (|12p to introduce the new problem: 

-^ + Jl Ul=P in [ °' Ll] - (14) 

It will replace (TTa]) in fJi£>. As evoked in Remark [TJ the reader can be easily convinced that 
it is particularly awkward to guess the value of L\. Indeed one has to specify the criteria that 
define 2-D effects, and in practical situations we may only be able to define L2 which is such 
that (p, n {x > L 2 }) C Q 2 d, or in other words L 2 > Ly. In this work we consider two different 
situations: 

- a funnel-shaped domain (see Figure ^ with a thin left part, so that we anticipate 2-D 
effects on the right (wide) part of the domain. In this case the definition of L 2 is based on 
a geometrical criterion. 

- a rectangular domain (see Figure [3} with a small aspect ratio e = H/ L (so that we can 
anticipate weak 2-D effects), but with some 2-D forcing terms occuring in the right end of 
the domain. In that case the definition of L 2 is based on the support of forcing terms. 

At this point we have defined an upper bound L 2 > L\, but the exact value of L\ remains 
unknown. From now on we choose an interface Lq without any a priori information (other than 
Lq < L 2 ) and decide to consider the model reduction (jT5J) on Hi — il PI {x < Lo}, while we keep 
the 2-D model in ^2 = fl {x > Lq}. Finally we have the two following systems 



d 2 u\ K 
' H 

w i(0) = 7i- 



I-D model: { ~ ^ Jj Ul - ^1 ^ (0, L ), ( 15 ) 
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Figure 2: Funnel-shaped computational domain. The domain is shallow for x < L%. 



Ii 




x = £a 



Figure 3: Rectangular computational domain 0. The domain is shallow: H/L <C 1, and we 
assume that the forcing terms are supported in {x > Inl- 



and 



2-D model: 



-Am 2 

du2 n t-2 
— — = on 1 rp 

OU 2 2 

—— + ku 2 = on Y B 

On 

U 2 



F 2 in fl 2 , 

r T ndn 2 , 
r B ndn 2 , 

72 on T_r. 



(16) 



where F\ = F and F 2 = -Fjn 2 . 
Two cases may occur: 

- Favourable case: Lq < L\, so that ill C and the following model reduction is relevant. 
In particular, hypothesis © holds so that Theorem [1] applies, 

- Unfavourable case: Lq > L\ : so that Oi <f_ fl\D and the 1-D model will not be able to 
reproduce the 2-D reality (in particular, hypothesis (O does not hold). 
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i = La 




Figure 4: Computational domains for the 1-D/2-D reduced model 



We now want to evaluate this model reduction in the favourable case Lq < L\, 

3 Coupling algorithm 

Let us consider the two models (|T5|) and (fll)]) to be coupled respectively through the interfaces 
x = Lq and T as shown in Figure 2) 

In coupling problems, the first difficulty lies in defining the coupling notion by itself, i.e. defining 
the quantities or values to be exchanged between the two models through the coupling interfaces. 
In our case and from a physical point of view, one may propose the following conditions, see [11] . 



which correspond to the conservation of u and its flux through the interface. 
Unfortunately these two constraints do not allow the well-posedness of the 2-D model and of the 
coupled problem. They are called defective boundary conditions in the literature, [2J, |14j . |12| . 
One should then rather apply a coupling method ensuring the following points: 

(i) the well-posedness of the 1-D and 2-D models 

(ii) the physical constraints are satisfied 

(iii) the control of the difference between the coupled solution and the reference one (corre- 
sponding to the 2-D model over the whole domain O). Indeed, due to the nature of the 
problem, the reader can be easily convinced that one does not expect to end with a solution 
of the 2-D model equal to the restriction of the reference one on f^- 

The coupling problem with (|17ap and (|17b|) has been studied in [TT] and [TJJ using variational 
and algebraic approaches. In this section, we propose an iterative coupling method based on 
classical Schwarz algorithms. These iterative methods were used for the first time in the context 
of dimensionally heterogeneous coupling in [5] and [5] to study a nonlinear hyperbolic coupling 




(17b) 



(17a) 



Inria 



DDM for a dim-heterogeneous problem 



9 



! 




Figure 5: Computational domains for the dimensionally homogeneous coupling problem 



problem. 

We will prove the convergence of these algorithms given an appropriate choice of boundary con- 
ditions at x = Lq and on T. Then we will study the solutions obtained after convergence and 
compare them to the global reference solution u defined by ([2]). Finally we will give some results 
regarding the choice of the coupling interface position. 



3.1 Schwarz algorithms 

Let us introduce first the Schwarz algorithms in the context of dimensionally homogeneous cou- 
pling. Consider the two systems, defined on fli and Q 2 shown in Figure [S] 



C u u = f u in SliCR" ( C v v = f v in n 2 C R" 

B° utu = 9u on dn i ut I B v uty = 9v on dVt° 2 ut 



The operators C u and C v are different. We assume that u and v have to satisfy the following 
constraints derived from the physics : 

du = C 2 v (19a) 
C[u = C' 2 v (19b) 

through the interface T, where Ci, C[, C 2 and C 2 are differential operators. 
To couple these two models we can implement the following iterative algorithm: 
For a given v° and at each iteration k > 0, solve: 

C u ii k+l = f u in Oi ( C v v k+1 = f v in VL 2 

B° ut u k+1 = g u on dn<[ ut then I B° ut v k+1 = g v on dQ^ 
B u u k+1 = B v v k on T [ B[,v k+1 = B' u u k+1 on T. 

Once convergence is achieved, the physical constraints (|19a[) and (|19b|) have to be satisfied. Then 
care should be taken to choose the operators B u , B' u , B v and B' v in order to ensure convergence 
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toward the unique solution denned by (fI5]l. (jTOa|) and (|T9bj) . see [J5] and [T5] . 

This method can be generalized to the dimensionally heterogeneous coupling case. Let assume 

that we have to solve the following 1-D model/2-D model coupled problem: 

=/i in fiiCR , [ £2^2 = f 2 in fl 2 C 

and 



B° ut Ul = g x on 30°"* \ B° ut u 2 = g 2 on dn° ut 

and we suppose that we have the following coupling constraints to satisfy at x = Lq and on T: 

C 1 u 1 (L ) = C 2 (TZu 2 )(L ) (20a) 
C' 2 u 2 (Lo,z) = C' 1 (£u 1 )(Lo,z) on T (20b) 

7lu 2 is a restriction of u 2 at x = Lo and £ui is an extension of iti (Lq) all along I\ More generally 
we can define the operators 1Z and £ as in [11] by: 

7e : A — ► A 
u 2 |r 1 — > Ku 2 \r 

£:A — > A 

u l|x=L 1 ^ ^ w l|x=L 

The spaces Ao et A are the trace spaces on the interface x = Lq for 1-D functions and on the 
interface T for 2-D functions. As mentioned in [TT], these two operators are not invertible. 
One may thus implement the following algorithm: 

For a given u 2 and at each iteration k > 0, solve: 

.fc+i — r r> ( r „,M-l 



and 



C 1 ul +L = fx in O x f £2^2 = /a in ^ 



Bl ut u k + l = 9l on 90? u * then < B^u^ 1 = g 2 on dVL° 2 ut 



Bxu^ 1 = B 2 Ku\ at i = L [ B' 2 u 2 +1 = B[£u k+1 on T. 

In practice, we do not have conditions such as (|20a|) (at 2; = Lq) and (j20b|) (for (x,z) on T), 
but only conditions at x = Lq such as (|17a[) and (|17b>[) . This leads us to make a choice of the 
operators 1Z, £, C±, C 2 , C[ and C 2 . 

In general the choice of the restriction operator 1Z is more straightforward than the choice 
of the extension operator £ one. In this study, since the 1-D model is obtained after some 
approximations and by averaging the 2-D model, it is reasonable to define 7Z as the vertical 
average. On the other hand, the question of the choice of the operator £ remains open. 
In [TT] and [T2], authors proposed a constant extension of (|17ap and (|17b[) along T, and impose 
the following strong coupling constraints: 

(1 f H 
Ul (L a ) = — u 2 {L 0l z)dz 

a if ^° 

ou 2 dm 



dx H J Q dx ' — — (L ,z) = ——(L ) on V 

v ox ox 

It is a choice among many others. One may also choose a multitude of operators Ci, C 2 , C[ and 
C 2 ensuring the physical constraints to be satisfied. 

The strategy that we adopt here is to choose, due to relations (TTQ)) and (fT2")) . a constant extension 
of iti along r and then to implement a family of Schwarz algorithms with appropriate boundary 
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conditions at x = Lq and on T. In this case Schwarz coupling algorithm reads: 
For a given u 2 and at each iteration k > 0, solve : 

d 2 u\ +1 



-u k+1 



and then solve 



dx 2 H 
^' +1 (0)=7i 
{ B 1 u k+1 {L a ) = B 1 u k 



-Au k+1 = F 2 in fl 2 



— = on 1 T 

on 



F x m (0,L ) 



(21) 



du k 2 +1 



dn 



ku\ +x =0 on r| 



(22) 



,fe+i 



= 72 on T R 



B 2 u k+1 = B 2 u k + l on T 

The linear operators B\ and B 2 will be defined such that the points (i), (ii) and (iii) (see 
introduction of Section are satisfied and such that the algorithm converges. 
We will first study the convergence of the coupling algorithm. Subsequently we move on the 
point (iii). 

To ensure the convergence of Schwarz algorithms in the case of classical domain decomposition 
without overlapping, it is proposed in |17| to use Robin operators. We will extend the use of 
these operators to our coupling problem. We define the operators B\ and B 2 for a given A > 
as follows: 



Bi 



d 
dni 



and 



B 2 = 



dn 2 



Xld 



Md 



(23) 



(24) 



where n\ and n 2 are the outward unit normal to the 1-D and 2-D domains respectively. We note 
that the operators B\ and B 2 ensure the well-posedness of the problem at each iteration. Let us 
study the convergence of Schwarz algorithm with this family of operators. 

3.2 Algorithm convergence 

Proposition 1. For each A > 0, the (u\ : u 2 ) Schwarz algorithm converges in iJ 1 (fii) x H 1 ^^ 
to (u*,u 2 ) that satisfies the physical constraints H7a)) , HI 7b)) . 

Proof: 

Let define the differences between two successive iterations: 



,fc+i i 



(x)=u k+1 (x)-u1(x), Vxe(0,L o ) 
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and 



e^ +1 (x, z) = v% +1 (x,z) - v%(x,z), V(x, z) G Q 2 



These functions satisfy the following systems: 
( 3 2 e k+1 



^e k+1 =0 in (0,L„) 



dx 2 H 
e k+1 (0) = 

de k ; +1 



dx 



(L ) + \e k+i (L ) 



^(L )+Xe k 2 (L ) 



and 



-Ae k+1 = in fl 2 



de k 2 +1 2 

— ^— = o on r| 

on 



Bp k+1 

° 2 ' K e k+1 = on r| 



9ti 



on r K 



de k+1 



-^2_(L ,z) + Xe k+1 (L ,z) = -^—(L Q ) + Xe k+1 (L Q ) on I\ 



The first two equations of (|2~5T) lead to: 

ej +1 (x) = a k+ i sinh(aa;), Va; G (0, L ) 

where otu+i G R and a = 



(25) 



(26) 



(27) 



If we take the vertical average of the boundary condition on V in (|26p , and due to the boundary 
condition at x = Lq, we obtain : 



3e k+ 



dx 



(L ) + \e k+1 (L ) = -^—( Lo ) + Xe k+1 (L ) 



f)pk+l q k+2 

dG2 -(L Q ) + Xe k+ \L ) = ^- iLo ) + \e k+2 (L Q ). 



dx 



This implies: 



and 



e k+1 (L ) 



2X 



(Aa k +2 + Ba k+ i) 



1 



de k+1 

q x (Lq) = | (Aa k+2 - Ba k+1 ) 



(28) 



(29) 
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where A = acosh(aLo) + Asinh(aLo) and B = — acosh(aLo) + Asinh(aLo). 

Now by multiplying the first equation of (|26p by e\ +1 and by integrating in 0,2, we obtain: 



\S7el +1 \ 2 dxdz 



on. 



-^—e^da = 
on 



then using the boundary conditions on T 2 ^ and T^, we obtain: 



We k+1 \ 2 dxdz + I kK + T&i 



dx 



(L Q ,z)e k+1 (L a ,z)dz. (30) 



We replace — ^ — (Loj-z) in fl30l) by its value obtained by using Robin boundary condition on V: 

ox 



\7e k+1 \ 2 dxdz + / K\e^ +1 \ 2 dx = 



del 



k+l 



(L ) + Ae* +1 (Lo)-Ae* +1 (L ,2) I e« +l {L ,z 



k+li 



,fc+l / 



(9.r 



Ba fc+1 F4- +1 (L )-A / |e* +1 | 2 (L ,z)dz 
Jo 

Ba k+1 H 



2A 



(Aa/c+2 + Ba i+ i) - A / \e k 2 +1 \ 2 {L , z)dz. 



We now replace e£ +1 (L , z) in ([3D) 

using the same Robin boundary condition: 



Ve k+i \ 2 dxdz + / K\e k+i \ 2 da = - 



1 / del 



k+l 



A 



dx 



-(L ) + Xe k+i (L ) 



— {L ,z)dz 



1 

A7o 



H 



(Ac 



dz 



A 

Ba k +iH 
2A 



<9x ^ ° A J 



H 



de k+1 



dz 



(Aa k +2 - Batk+i) - -r 



H 



de k+1 



dz. 



Due to the fact that A > 0, we deduce that: 

Ba k+1 H 



and 



Thus: 



2A 

Ba k+ iH 
2A 

l2„2 



{Aa k+2 +Ba k+1 ) > 
(Aa k+2 - Ba k+1 ) > 0. 



A'a 2 k+2 - B 2 a 2 k+1 < 



Then we obtain: 



l k+2 



< 



B 2 



— acosh(aio) + Asinh(aLg) 



acosh(aio) + Asinh(aLo) 



< 1 
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< 



and finally 

Ctk+2 
"fe+1 

So that the sequence (ctk)keft converge to zero. 

Let us now remark that for all k > 0, n > 0, we have: 



< 1. (32) 



,fc+n 



n-1 

J2 e1 +p+1 

p=0 



and 



d (u k+n - u k ) ^ de\ +p+1 

dx dx 

p=0 



Using the relation (p?T)) and the fact that the sequence (ak)k^n converges, we can prove that 
du\ 
dx 



du k 

(ui)fe 6 N and (— — )/cgN are Cauchy sequences in L 2 (ili). So that (iti)fceN is a Cauchy sequence 



in ff^fix). 

In the same way we observe that for all fc > 0, n > 0, we have: 

n-1 

V(u k 2 +n -u k 2 ) = ^Ve^ +1 

p=0 

Using (|3"Tj) . we deduce that: 

/ \\7e k 2 +1 \ 2 dxdz < Bak + lH ( A a k+2 + Ba k+1 ) 

and then we can prove that (Ve^fceN is a Cauchy sequence in L 2 (ft 2 ) and due to the Poincare 
inequality we have also (tif)fceN is a Cauchy sequence in L 2 (il 2 ). So that (lif)fceN is a Cauchy 
sequence in H 1 (£l 2 ). 

To conclude we have prove that (u k ,u 2 ) Schwarz algorithm converges in x H 1 ^^). 

Moreover, at convergence the limit (ui,u 2 ) verifies B 2 u 2 = B 2 u\ and B\u\ = Biu 2 . Taking the 
vertical average on T gives two linear combinations of the constraints (|17aj) . (|17bp .D 

Remark 2. 



• The Schwarz algorithms converge for all A positive, but we remark that for A = a coth(aLo)> 
we have exact convergence in two iterations. Indeed we have in this case: 

de k+1 

dx~ + acot H aL o) e i +1 = Vfc > 0, 

and then: 

de k+1 de k+1 

B 2 e k+1 = ^— + acoth(aL )e k+1 = ^— + a coth(aL )ef +1 =0 Vfc > 0, 

dx dx 

d 

The operator — \- a coth( a Lq) Id corresponds to the absorbing operator of the 1-D model. 

dn\ 

We denote by X opt this value of A, see for example JJEj. 
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If we take Ai ^ X2, we have a priori (u* 1 , W2 1 ) 7^ (u^ 2 , it^ 2 )- This is in accordance with the 
ill-posedness of coupling problem defined by H5\) . \lb}) . H7a\) and \17b\) . 



For the sake of clarity we will denote (ux,U2) the limit of Schwarz algorithm instead of (u^u^) 



3.3 Control of the difference between the coupled solution and the 
global reference solution 

Unlike the case of domain decomposition, at convergence of Schwarz algorithm, we have it 2 ^ it |n 2 
due to the model reduction. But as mentioned above, we have chosen the family of Robin 
operators in order to get some control of the difference between U2 and u\q 2 . In fact we have the 
following result: 



Theorem 1. for each A > 0, let (111,112) denotes the limit of the Schwarz algorithm. If Lq < L\ 
then there exists M(A) > such that: 



lhn 2 -Msnch) < M(X)e\/l + 5 2 



(33) 



where S = 



Li 



L\ — Lq 



Proof 

The function u\n 2 — U2 is the solution of the system: 



— A(w — U2) = in O2 

d(u - u 2 ) 
dn 

d(u - u 2 ) 



on H 



dn 

u — \i2 = on 



niu — U2) = on T\ 



By multiplying the first equation by u — U2 and using the boundary conditions on Tj, U T 2 B UTr, 
we obtain: 



/ V (u — U2) \ 2 dxdz + / k\u — U2\ 2 dx — I 
J no Jr 2 „ Jr 



d(u - u 2 ) 
dn 



(u — U2) dz = 0. 



(*) 
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du du 

The integral term on T is reformulated using the boundary condition — \-\1i2 = — r h Aui 

dx ox 

satisfied by the limit u 2 : 

d{u~u 2 ) . . f H d(u-u 2 ) 

(u~u 2 )dz = (Lo,z) (u - u 2 ) (L ,z)dz 

On J Ox 

rB du 



(L , z) (u~u 2 ) (L , z)dz 
H ( du 

f ~~g^~( L °) + ^ u i(Lo) - Xu 2 (L , z) ) (u - u 2 ) (L , z)dz 
1 ' du 

~g^( L o, z) + Xu(L , z) ) (u - u 2 ) (L , z)dz 
H ( du 

~~jfa( L °) + ^"lOM ) {u - u 2 ) (L , z)dz 
-X I (it — u 2 ) (L , z)dz. 



The first term reads: 

H / du \ , f H ( du 



-g^+Xuj (Lo,z)(u-u 2 )(L ,z)dz = / I -— + Ait ) {L ,z) (u(L , z) - u{L Q )) dz 



H 







H 



du 

-— + Xu ) (L , z) (u(L ) - Ui(L )) dz 
du 

-— + Xu ) (L , z) (ui(L ) - u 2 (L , z)) dz. 
du 

Due to the relations ([TU| and (fT2"jl and to the fact that — — (L ,z) + Xu(L Q ,z) = 0(1), we 
deduce: 

( - — + Xu) (L , z) (u(L , z) - u(L )) dz = 0(e 2 ) 

In the same way, if we assume that Lq < Li, so that 2-D effects are insignificant in Six n 
{Lq < x < Li}, and applying a similar asymptotic analysis as in the first section to the 2-D 
model defined in VL 2 , we can deduce that: 

u 2 (L Q ,z) = u 2 (L )+O(8 2 e 2 ) 

= Ul (L )+O(S 2 e 2 ), Vze[0,ff] (34) 



So that: 



Finally: 



du 



— + Xu ) (L ,z) (wi(L ) - u 2 (L , z)) dz = <3(<5V) 



(~£ + Xu ) (L °' z) (s(io) ~ ui(Lo)) dz=H + xu 
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• Since ui(Lq) = u 2 (Lq), the second term reads: 

^(Lo) + A«i(io)) (u - u 2 ) (L , z)dz = -H (- ^"g"^ + A(«i - u)) (L„) (« - «i) (L ) 

-h(-P+\u)(L )(u- Ui )(L ) (**) 



9x 

We reformulate the first term of the right. Note that the function u\ — u satisfies the equation: 
d 2 (u± - u) 



dx 2 



-(x) + a [u\ — u)(x) = a (u(x) — u(x, 0)), Vx G (0, Lq) 



So that by multiplying this equation by u\ — u, after integration on (0, Lq) and use of the 
boundary condition ui(0) = u(0), we obtain: 

J ° f d(Ul ~ ^\x)dx + a 2 f L ° ( Ul - S) 2 - ^V^ ol^i - u)(L ) = 



dx J Jn dx 

r-Lo 

2 



a (u(x) — u(x, 0))(ui(x) — u){x)dx 



thus: 



d{u\ — u) f L " 

(Lo)(ui — u)(Lq) = / a 2 (u(x) — u(x, 0))(ui(x) — u(x))dx — A\{u\ — u, u\ — u) 

ox J 

f L ° fd(u\ — u)\ 2 f L ° 2 

where A\ (u\ — u, U\ — u) = I ( ) (x)dx + a 2 / (ui — u) (x)dx. 

Jo \ ox J J 

And then becomes: 

' dui <- L 
dx 

—HA\{ui — u, ui — u) + XH (u — ui) 2 (Lq) 
-h(-^ + Xu](Lq)(u-ui)(Lq) (36) 



— — (Lq) + \ui(Lq) I (u — U2) dz = Hi a 2 (u(x) — u(x, 0))(ui(x) — u(x))dx 



• To recap, the boundary term on T in Q becomes: 



U2) (u - u 2 ) dz = 0(e 2 ) + 0(S 2 e 2 ) + H [ ° a 2 (u(x) - u(x, 0j)(ui(x) - u)dx 
on J 

—HA(ui — U, Ui — u) + XH (u — ui) 2 (Lq) 

—A / (u — u 2 ) dz 
Jo 

We first observe that: 

Lq a 2 rL a 2 rL 

a 2 (u(x) — u(x,0))(ui(x)—u(x))dx<— / (u(x) — u(x , 0)) 2 dx + — / (ui(x) — u(x)) 2 dx. 
2 J 2 Jo 
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It follows that: 

,2 r L o 



I d{u_U2)_ u ^ C {l + S 2 )e 2 ) + H— [ \u(x)-u(x,0)) 2 dx 

Jr on 2 J 



2 f-L 



—HAi(ui—u,ui—u) + H— J (ui(x) — u(x)) dx 



o 

\2 l t \ \ \ I \2 



+XH(u — ui) (L ) - A / (u — u 2 ) (L ,z)dz 



ii 



where C is a positive constant. 
Then we have: 



a 2 r La 



~A\{ui~u,ui—u) + — j (ui(x) — u{x)) 2 dx < 



and finally, using the definition of u(Lq) and the relation u\(Lq) = u 2 (Lq), we obtain: 

\ 2 













H 








{I 










a! 


U 


H 





ldz / (u — u 2 ) 2 dz — A / (it — U2) 



< 0. 

We now come back to which gives: 

|V (u - u 2 ) \ 2 dxdz + / n\u-u 2 \ 2 dx < M(l + S 2 )e 2 
and thus: 

ST2\_2 



/ |V(u-u 2 ) | 2 da:dz < M{l + 5 2 )e 2 



Where M denotes a positive constant. 

Finally, due to the fact that u — u 2 = on r#, and by using Poincare inequality we can deduce 
the inequality (|55|) .D 



Remark 3. 



This proposition fails if we choose the interface position in a zone where 2-D effects are 
significant. In this case relation {3$ is no more available. 

The right term of k33]) is also an upper bound of ju^ 1 — w^ 2 f or °^ an ^ ^ 2 

positive. 
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4 Numerical results 

The test cases presented in this section illustrate the coupling method of 1-D and 2-D elliptic 
equations based on Schwarz algorithm. All the computations have been done using the software 
package Freefem++ [TH], with a P2 finite element discretization. 

In the first part of this section, the two test cases will be described in details. In the second and 
third parts, we will focus on one hand on the Scharwz algorithm convergence and on the other 
hand on the comparison of the coupled solution with the reference solution in order to enlight 
the theoretical results obtained in the previous paragraphs. 

4.1 Description of the test cases 
4.1.1 Test #1: 

The first test case is concerned with the solution of the 2-D problem (J2J) where the domain is 
a rectangle ft = [0,L] x [0,H] which is assumed to be uniformly shallow: H <C L. Let us 
consider that the right-hand side term F(x, z) of the full 2-D problem is F(x, z) = mexp(— (x — 
2ttz 

x*) 2 ) sin(— — - ), where x* < L. 
H 

The global reference solution u is displayed in Figure [51 




-0.0949 -0.000242 0.0944 



Figure 6: 2-D reference solution for the first test case where L = 20, x* = 19, H = 0.5 and 
k = 0.001. 



We notice that the 2-D effects are due to the particular form of the forcing term F and are 
located around x* . 

Now let us define the coupled model. 
The interface is located at x = Lq < x* as shown in Figure [7J In the part of the domain 
ill = [0, Lq] x [0, H], we assume a priori that the 2-D effects are negligible and consequently we 
replace the full 2-D equations by the 1-D model (see (TT5"}) ). 



4.1.2 Test #2: 

In this second test case, the 2-D effects are due to the funnel-shaped geometry of the domain 
(see Figure I5a|) . and the forcing term is constant (= 1). 

The reference solution in the whole domain is displayed in Figure |H1 
The coupled model is defined by splitting the domain in two parts. The interface T is located at 
x = Lq, < Lq < L>2 as shown in Figure |S] 
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.1 ■ 



L 



(a) Computational domain for the 2-D reference model 



V. 



(b) Computational domain for the 1-D/2-D reduced model 



Figure 7: Computational domains for both the reference and reduced models in test case $1. 
For the reduced model (b), the 1-D/2-D interface T is located in x = Lq. 



4.2 Convergence of the Schwarz algorithm 

In this section we provide numerical results to assess the theoretical results of JJl We are inter- 
ested in illustrating the optimal convergence of Schwarz algorithm for the parameter A = X op t- 
Figure HOal shows the difference between the iterates of the Schwarz algorithm in L°° norm for 
the two test cases. 

As demonstrated in JJJ the Schwarz algorithm converges in two iterations for the optimal parame- 
trer A = X op t ■ It is important to notice that this result is independent of the interface location. 



4.3 Difference between the coupled solution and the full 2-D solution 

One important point in the analysis of the accuracy of the coupling procedure is the comparison 
of the coupled solution with the reference solution as a function of the interface location x = Lq. 
Contrarily to the classical domain decomposition problems, here there is a difference between 
the (converged) coupled solution and the full 2-D solution; this difference is due to the model 
reduction that is performed in the I-D part of the domain. This difference depends on the loca- 
tion chosen to discriminate between 1-D and 2-D regions. Figures fllal and f l lbl left show the H 1 
error between coupled and reference solutions as a function of the interface location for the two 
test cases. Figures lllal and If f bl right show the H 1 error in fl2 between coupled and reference 

solutions as a function of e — — for the two test cases. 

It is interesting to notice that for both test cases there is a discontinuity in the curve rep- 
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l + H 



L2 



(a) Computational domain for the 2-D reference model 




(b) Computational domain for the 1-D/2-D reduced model 



Figure 8: Computational domains for both the reference and reduced models in test case #2. 
For the reduced model (b), the 1-D/2-D interface V is located in x = Lq. 



resenting the error as a function of Lq (see Figure 111! left column) . This discontinuity occurs 
both for the numerical difference between U2 and uq 2 (black curve), and the theoretical curve 
(in red) corresponding to the right-hand-side of estimate (|33p . Indeed, if Lq is greater than a 
certain threshold, the error grows very rapidly (and S — > 00 in estimate (|33|) ). This could be 
an indication of the real (a priori unknown) value of L\ (see discussion at the end of Section [5] 
above) . 
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Figure 9: 2-D reference solution for the second test case where L = 2, H = 0.05, I = 3 and 
k = 0.001 




1 = 3 and re = 0.001. 

5 Conclusion 

In this paper we studied a linear boundary valued problem set in a 2-D domain, and assume 
that the solution may be approximated by a 1-D function in some part of the computational 
domain. We thus derive a reduced model that consists coupling a 1-D model (wherever we think 
it is legitimate) together with the original 2-D system (everywhere else). 

The model reduction is performed thanks to a small aspect ratio hypothesis, with an integration 
in the shallow direction (we mimic the derivation of the shallow water equations). After this 
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(b) Test case # 2: (left): L = 2, H = 0.05 and re = 0.001. (right): L = 2, L = 1.5 and re = 0.001. 



Figure 11: Relative error as a function of Lq (left) and e (right) between the coupled solution 
and the 2-D reference solution in test case #1 (top) and #2 (bottom). In the left column, the 
red curves correspond (for both test cases) to the RHS of estimate (f3"31 . 



derivation we introduce an iterative method that couples the 1-D and 2-D systems and we prove 
some convergence results. One original aspect of this work is the particular attention that is paid 
to the location of the 1-D/2-D interface. These theoretical results are illustrated with numerical 
simulations that underline the importance of the interface position, but also the way 1-D and 2-D 
models are coupled (boundary conditions at this interface). All these aspects, that have been 
studied here with a linear model, will be considered in a forthcoming study of dimensionally 
heterogeneous modelling in fluid dynamics. 
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